Breakdown of Landau-Ginzburg- Wilson theory for certain quantum phase transitions 
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5-1 '. 

Oh. The quantum ferromagnetic transition of itinerant electrons is considered. It is shown that 

' the Landau-Ginzburg- Wilson theory described by Hertz and others breaks down due to a singular 

coupling between fluctuations of the conserved order parameter. This coupling induces an effective 
long-range interaction between the spins of the form l/r^"*"^. It leads to unusual scaling behavior 
at the quantum critical point in 1 < d < 3 dimensions, which is determined exactly. 

"r-| '! PACS numbers: 64.60.Ak , 75.10.Jm , 75.40. Cx , 75.40. Gb 

o 

1^ ' One of the most obvious examples of a quantum phase transition is the ferromagnetic transition of itinerant electrons 
at zero temperature T as a function of the exchange coupling between the electron spins. Hertz [Q derived a Landau- 
Ginzburg- Wilson (LGW) functional for this case in analogy to Wilson's treatment of classical phase transitions, and 

■ analyzed it by means of renormalization group methods. He found that the critical behavior in dimensions d = 3, 2 
c/2 is mean-field like, since the dynamical critical exponent z decreases the upper critical dimension compared to the 

classical case. In a quest for nontrivial critical behavior. Hertz studied a model with a magnetization confined to c? < 3 
dimensions, while the coefficients in the LGW functional are those of a 3-d Fermi gas. For this model he concluded 
that d^ — 1, and performed a 1 — e expansion to calculate critical exponents in d < 1. Despite the artificial nature of 
his model, there is a general belief that the qualitative features of Hertz's analysis, in particular the fact that there is 
(~| mean- field like critical behavior for all d > 1, apply to real itinerant quantum ferromagnets as well. 
O In this Letter we show that this belief is mistaken, since the LGW approach breaks down due to the presence of soft 

1^1 ^ modes in addition to the order parameter fluctuations, viz. spin-triplet particle-hole excitations that are integrated 
out in the derivation of the LGW functional. These soft modes lead to singular vertices in the LGW functional, 
invalidating the LGW philosophy of deriving an effective local field theory in terms of the order parameter only 
In Hertz's original model this does not change the critical behavior in d > 1, but it invalidates his 1 — e expansion. 
More importantly, in a more realistic model the same effect leads to nontrivial critical behavior for 1 < d < 3, which 
we determine exactly. 

■ Our results for realistic quantum magnets can be summarized as follows. The magnetization, m, at T = in a 
T-H ' magnetic field H is given by the equation of state 

^ . 

tm + vm'^ + um^ = H, (1) 

^ : where t is the dimensionless distance from the critical point, and u and v are finite numbers. From (Q) one obtains 
' the critical exponents /? and S, defined by m ~ and m ~ H^^^ , respectively, at T = 0. For /3 and 6, and for the 
. correlation length exponent i^, the order parameter susceptibility exponent r/, and the dynamical exponent z, we find 

c ■ 

Q, f3^iy=l/{d~l), ri^3-d, S = z = d, (1< d < 3), (2) 

O . 

' and j3 = V = 1/2, ri = Q,6~z^'i for d > 3. These exponents 'lock into' mean-field values at d = 3, but have 

• ^ , nontrivial values for d < 3. In d = 3, there are logarithmic corrections to power-law scaling. Eq. (|l|) applies to T = 0. 

^ ' At finite temperature, we find homogeneity laws for m, and for the magnetic susceptibility, Xm, 

. . m(t,r,i/) =6-^/''m(^6l/^^&*/^i^6*^/'') , (3a) 

Xm{t,T,H)^}p/-Xm{th^''',Tb^''',Hb'^'n , (3b) 

where b is an arbitrary scale factor. The exponent 7, defined by Xm ~ at T = H = and the crossover exponent 
(j) that describes the crossover from the quantum to the classical Heisenberg fixed point (FP) are given by 

7 = /3(<5-l) = l , (^ = 1^ , (4) 
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for all d > 1. Notice that the temperature dependence of the magnetization is not given by the dynamical exponent. 
However, z controls the temperature dependence of the specific heat coefficient, 7y = cy/T, which has a scale 
dimension of zero for all d, and logarithmic corrections to scaling for all d < 3 

jv{t,T,H)^e{3-d)\nb + -fvitb^^'',Tb'',Hb^'^^'') . (5) 

Eqs. (|l|) - (^) represent the exact critical behavior of itinerant quantum Heisenberg ferromagnets for all d > 1 with 
the exception of d = 3, where additional logarithmic corrections to scaling appear. We are able to obtain the critical 
behavior exactly, yet it is not mean-field like. The exactness is due to the fact that we work above the upper critical 
dimension d^ = 1. The nontrivial exponents are due to a singular coupling between the critical modes which leads, 
e.g., to the unusual term ^ v m (^. Experimentally, we predict that for 3-d magnets with a very low Tc there is 
a crossover from essentially mean-field quantum behavior to classical Heisenberg behavior. In d = 2, where there 
is no classical transition, we predict that with decreasing T, long-range order will develop, and the quantum phase 
transition at T = will display the nontrivial critical behavior shown above. 

We now sketch the derivation of these results. A more complete account of the technical details will be given 
elsewhere [Q. We consider a d-dimensional continuum model of interacting electrons, and pay particular attention to 
the particle- hole spin-triplet contribution to the interaction term in the action, Sfu^, whose (repulsive) coupling 
constant we denote by J. Writing only the latter explicitly, and denoting the spin density by ns, the action reads, 

S = So + ^So + (J/2) J dx n,{x) ■ n,(a;) , (6) 

where Sq contains all contributions to the action other than Sf^^^. In particular, it contains the particle-hole spin- 
singlet and particle-particle interactions, which will be important for what follows. J dx = J dx /q^^"^ dr, and we use 
a 4- vector notation x = (x, t), with x a vector in real space, and r imaginary time. Following Hertz, we perform 
a Hubbard-Stratonovich decoupling of 5**^^^ by introducing a classical vector field M(a:) with components that 
couples to ns(a;) and whose average is proportional to the magnetization, and wc integrate out all fcrmionic degrees 
of freedom. We obtain the partition function Z in the form 

Z = e-''"/'^ J D[M] exp[^^M]] , (7a) 

where Fq is the noncritical part of the free energy. The Landau-Ginzburg- Wilson (LGW) functional <I> reads 

$[M] = i /" dxdy j6{x-y)M{x)-M{y) + f2an f dxi . . . dxn x^^LM, . . . ,x,,)M'' {x^) . . . AT- {xn) , (7b) 



where a„ — (— l)"+^/n!. The coefficients x^"-* in (''b) are connected n-point spin density correlation functions of a 
reference system with action 5*0 . The particle-hole spin-triplet interaction J is missing in the bare reference system, 
but a nonzero J is generated perturbatively by the particle-particle interaction contained in Sq . The reference system 
then has all of the characteristics of the full action S, except that it must not undergo a phase transition lest the 
separation of modes that is implicit in our singling out Sf^^ for the decoupling procedure breaks down. 

X^^^ is the spin susceptibility of the reference system. Performing a Fourier transform from x = (x, t) to q = (q, $7) 
with wavevector q and Matsubara frequency Vt, we have for small q and Q. |^, 

x(2)(q,l])=Xo(q)[l-|f^l/|q|] , (8a) 

where q and O are being measured in suitable units, and Xo(q) is the static spin susceptibility of the reference system. 
We now use the fact that in a Fermi liquid at T = 0, xo is a nonanalytic function of q of the form 

Xo(q->0)- const -Iql'^-i-q^ . (8b) 

Here we have omitted all prefactors, since they are irrelevant for our purposes. This holds for 1 < d < 3; in d = 3 the 
nonanalyticity is of the form q^ In |q| 0. Using (||), and with = T X^jjj: tJ^c Gaussian part of $ can be written, 

$(2)[M]= /"M(q)[io + c„|q|'^-i+c,q2+Crf|0|/|q|]M(-g) . (9) 

Jq 

Here = 1 — TtX^^Hq ~* ^i^n — 0) is the bare distance from the critical point, and c„, Ca and Cd are constants. 
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For the same physical reasons for which the nonanalyticity occurs in ( |8b| ) , the coefficients x*^"^ in (^) are in general 
not finite at zero frequencies and wavenumbers. Let us focus in x'^^K which will be the most interesting one for our 
purposes. Again, standard perturbation theory shows that it is given schematically by [|| 

X^*) - const + v [ [|k| + |w„|] - M + v\pf-^ . (10) 

Jk 

Here we have cut off the singularity by means of a wavenumbcr |p|, and u and v are finite numbers. More generally, 
the coefficient of |M|" in $ for |p| behaves like x''"'' = u'") |p|'^+^^". This implies that $ contains a nonanalyticity 
which in our expansion takes the form of a power series in |Mp/|pp. 

The functional $ can be analyzed by using standard techniques We look for a FP where Cd and either c„ (for 
1 < d < 3), or Ca (for d > 3) are not renormalized. This fixes the critical exponents ry and z. Choosing the scale 
dimension of a length L to be [L] — —1, standard power counting j|] then yields the scale dimension of t''"-' to be 
= — (n — 2){d — l)/2. All non-Gaussian terms are thus irrelevant for d > 1, and they all become marginal in 
d = 1 and relevant for d < I. Several features of the critical behavior follow immediately. The critical exponents 77 and 
z are fixed by the choice of our FP, and v and 7 as given in (||) and (^ are obtained by considering the q-dependence 
of the Gaussian vertex (^. We determine the equation of state by taking the term of order |M|^ in $ into account. 
X^'*^ is dangerously irrelevant with respect to the magnetization. We have shown jj] that for scaling purposes the 
cutoff IpI in ( p^ can be replaced by m. From this and (9) we obtain the effective equation of state as given in (|^). 

These results completely specify the critical behavior at T = 0. Their most interesting aspect is the nontrivial 
exponent values found for 1 < d < 3, which can nevertheless be determined exactly. The reason for this is the 
|q|''~^-term in the Gaussian action (^. It reflects the fact that in an interacting electron system, static correlations 
between spins do not fall off exponentially with distance, but only algebraically like 7--(2<i-i), This slow decay leads 
to a long-range interaction in the effective action which falls off like l/r^'*"^, see (ph. The critical behavior of classical 
Heisenberg magnets with such a long-range interaction has been studied before 0]. 



We now turn to the T-dependence of the specific heat, cy. We expand the free energy functional (7b) about the 
expectation value, m, of M to second order, and then perform the Gaussian integral to obtain the partition function. 
The free energy is obtained as the sum of a mean- field contribution given by <i>[m], and a fluctuation contribution 
given by the Gaussian integral. The latter yields the leading nonanalytic term in the free energy. We find |Q] that 
effectively H and T have the same scale dimension, viz. d(= z), and that a,tt — Q there is a logarithmic T-dependence 
of 7y for alll < d <2). If we put the i-dependence back in, we obtain that the scale dependence of 7y is given by (|^). 

For the magnetization the leading T-dependence is given by the mean-field contribution to the free energy. We 
calculate the temperature corrections to the equation of state (|l|) and find that for m >> T (in suitable units) m'^ in 
(^ will be replaced by m'^[l -I- const x T/m -I- . . .], while for m << T, t is replaced hy t + T^/" . The effective scale 
dimension of T in m is therefore 1 [not z), and we obtain for m and Xm the homogeneity laws given by (|). Thus, the 
relevant operator T in (H) reflects the crossover from the quantum to the classical FP rather than dynamical scaling. 
Accordingly, we have written the T-dependence in (|^) in terms of a crossover exponent (p which is given by ^ . 

Next we briefly discuss Hertz's original model, which differs from the one discussed above in two ways. First, the 
reference ensemble consists of noninteracting electrons. Second, the coefficients x*^"-* are taken to be the correlation 
functions of a i-d fermion system, x^^^ is then simply the Lindhard function, so (^ gets replaced by, 

x(^Hq,f^) = i-q'-|^^l/|q| + --- • (n) 

Due to the missing interaction in the reference ensenble, x^'^\''\^ 0) is now analytic at |q[ = 0. The resulting quadratic 
term in (^) allows for a Gaussian FP with mean-field static exponents and a dynamical exponent z = 3 ||l[] . Whether 
this FP is stable depends on the higher x'""*- Hertz considered only the limit q = 17 = 0, where all of these terms are 
finite numbers and irrelevant for d > 1. The quartic term is marginal in d = 1 and relevant for d < 1 j^]. 

The striking difference between the finite coefficients in Hertz's model and the diverging ones in the realistic model 
above is due to the latter containing interactions in the reference ensemble. The interactions lead to frequency mixing, 
and hence to soft particle-hole excitations contributing to the x^"^ even in the limit of zero external frequency. A 
similar effect is achieved for noninteracting electrons by considering correlation functions at nonvanishing external 
frequency. Therefore we include the higher order terms in an expansion of the x^"-* in powers of and analyze the 
arising LGW functional by the same power counting arguments as above. The details of this calculation will be 
presented elsewhere [4] . We find that all non-Gaussian terms are still irrelevant for d > 2 and the critical behavior is 
mean- field like. In d = 1, however, the x*^""* change their functional form so that an infinite number of operators is 
relevant [not marginal) with respect to the Gaussian FP in d = 1 and below. Therefore the upper critical dimension 
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is not one, but rather the 1-d sytem is below its upper critical dimension, and will show critical behavior that is 
substantially different from mean-field behavior. 

We conclude with a few remarks. First, the vertices in the LGW functional discussed here are singular only if the 
order parameter is conserved, and only at zero temperature, which means that the phenomenon is confined to the 
quantum magnetic transition. Second, our conclusion, although derived for the special example of itinerant quantum 
ferromagnetism, is rather general: We expect the LGW formalism to break down whenever there are soft modes other 
than the critical order parameter fluctuations that couple to the order parameter. The general rule is that all of the 
soft modes should be retained on equal footing in the effective theory. If any of them are integrated out, the resulting 
penalty are ill-behaved coefficients in the LGW functional. This has been shown recently for disordered electrons . 
The present results indicate that the underlying principle is very general. Indeed, it also applies to classical phase 
transitions with additional soft modes. However, there are many modes that are soft at T = but acquire a mass at 
finite temperature, making quantum phase transitions more likely candidates. Finally, we mention that Sachdev pl| ] 
has noted that something must be wrong with Hertz's theory in d < 1, since it violates an exact exponent equality 
for quantum phase transitions with conserved order parameters. He suspected Hertz's omission of the cubic term in 
the LGW functional to be at fault. Our analysis provides instead the explanation given above, namely the presence 
of infinitely many relevant operators due to the soft particle-hole excitations. 

This work was supported by the NSF under grant Nos. DMR-92-17496 and DMR-95-10185, by the DAAD, by the 
DFG under grant number Vo 659/1-1, and by the NATO under grant No. CRG-941250. DB would hke to thank 
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